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ABSTRACT 


The paper presents a correction to the scaling relations for red-giant stars using model-based masses 
and radii. We measure radial-mode frequencies from Kepler observations for 3,642 solar-like oscillators 
on the red-giant branch and use them to characterise the stars with the grid-based modeling. We 
determine fundamental stellar parameters with good precision: the typical uncertainty is 4.5% for 
mass, 16% for age, 0.006 dex for surface gravity, and 1.7% for radius. We also achieve good accuracy 
for estimated masses and radii, based on comparing with those determined for eclipsing binaries. 
We find a systematic offset of ~ 15% in mass and ~ 7% in radius between the modeling solutions 
and the scaling relations. Further investigation indicates that these offsets are mainly caused by a 
systematic bias in the Av scaling relation: the original scaling relation underestimates Av value by 
~ 4%, on average, and it is important to correct for the surface term in the calibration. We find no 
significant offset in the Vmax scaling relation, although a clear metallicity dependence is seen and we 
suggest including a metallicity term in the formulae. Lastly, we calibrate new scaling relations for 
red-giant stars based on observed global seismic parameters, spectroscopic effective temperatures and 
metallicities, and modeling-inferred masses and radii. 


Keywords: star: oscillation 


1. INTRODUCTION 


Asteroseismology of red giants has been undergoing extensive development in the last decade, thanks to the space 
missions like CoRoT (Baglin et al. 2002), Kepler(Borucki et al. 2009), K2 (Howell et al. 2014) and TESS (Ricker et al. 
2015). High-quality photometry data for thousands of red giants are now available for asteroseismology. Automated 
data reduction tools for measuring global seismic parameters (e.g. SYD Pipeline; Huber et al. 2009a), acoustic modes 
(e.g. Pbjam; Nielsen et al. 2021), and mixed modes (e.g. FAMED; Corsaro et al. 2020) make it possible to study 
solar-like oscillations in large samples of stars. Compared with studying individual stars, a population study covering a 
wide parameter range has great advantages. For instance, the systematic characterisation of solar-like oscillations and 
granulation can constrain stellar masses and radii with good precision (e.g. Yu et al. 2018a). In addition, measuring 
gravity-mode period spacings for red giants allows us to precisely monitor their stellar evolution from the subgiant 
stage to the end of helium burning phase (Mosser et al. 2014; Vrard et al. 2016). Finally, core rotation periods, 
measured from the splitting in oscillation modes (e.g. Mosser et al. 2012a), show a significant spinning down during 
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the red-giant phase and provide constraints on the angular momentum transport between the core and the envelope 
(Aerts et al. 2019; Tayar et al. 2019). 

Studying a large sample of stars on the red-giant branch (RGB) could advance the resolution of puzzles at this 
evolutionary stage. One of these is the remarkable accuracy of the scaling relations. Scaling relations were introduced 
to predict the oscillation frequencies of solar-type stars (Ulrich 1986; Brown et al. 1991; Kjeldsen & Bedding 1995): 


Av (p 0.5 T 
Avo \po 
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~ , (2) 
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Here, Vmax is the frequency at the maximum amplitude of stellar oscillations, and Av is the mean frequency separation 
between consecutive modes of the same spherical degree. These relations can be inverted to estimate masses and radii 
(Stello et al. 2008; Kallinger et al. 2010a): 
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The solar reference values have been measured in previous studies with slight differences, depending on the data and 
the reduction method. In this work, we adopt Vmax,o = 3090+30uHz, Avo = 135.1+0.14Hz, and Tefo = 5777K 
(Huber et al. 2011). 

The accuracy of scaling relations has been widely discussed (Stello et al. 2009; White et al. 2011; Huber et al. 2012; 
Mosser et al. 2013; Li et al. 2021) and noticeable systematic offsets have been seen in RGB stars. Stars in eclipsing 
binary (EB) systems can provide direct tests of the scaling relations. For instance, Gaulme et al. (2016) studied ten 
high-luminosity red-giant stars in eclipsing binaries by comparing masses and radii obtained from the scaling relations 
against spectroscopic orbital solutions using dynamical modeling. They found that the scaling relations overestimate 
the radii by about 5% and the masses by about 15%. With the similar approaches, Brogaard et al. (2018) and Themeßl 
et al. (2018) also found overestimated masses and radii derived by the scaling relations. The radius measured based 
on interferometry provides another direct test. For instance, Huber et al. (2012) combined interferometric angular 
diameters, Hipparcos parallaxes, asteroseismic densities, bolometric fluxes and high-resolution spectroscopy to derive 
near model-independent fundamental properties for one subgiant and four red giant stars near the base of RGB. Unlike 
the EB results, their measurements showed good agreement with the scaling relations. 

Determining masses and radii from the detailed modeling is another way to verify the scaling relations. Pérez 
Hernández et al. (2016) characterised 19 low-luminosity Kepler red giants based on a minimisation procedure combining 
the frequencies of the p-modes, the period spacing of gravity modes, and the spectroscopic data. They found no obvious 
differences between inferred masses and radii from the modeling and the scaling relations. It should be noted that 
the apparent conflict between these various studies may relate to the evolutionary stage rather than differences in the 
methods. As presented by Sharma et al. (2016), the systematic offset of the Av scaling relation varies significantly for 
different evolutionary stages. Global diagnosis of the scaling relations that covers a wide parameter range gives a big 
picture of the systematic uncertainty in the scaling relations. Huber et al. (2017), Zinn et al. (2019), and Hall et al. 
(2019) used Gaia parallaxes to empirically demonstrate that evolution-dependent corrections to the scaling relations 
are necessary to improve the accuracy of asteroseismic radii for thousands of Kepler RGB stars. 

Various methods have been implemented to correct the scaling relations. A simple approach is changing the solar 
reference values of Ave and Vmax,o (Kallinger et al. 2010b; Chaplin et al. 2011; Guggenberger et al. 2016; Themefil 
et al. 2018), however, the method does not consider the change for different evolutionary phases. Sharma et al. 
(2016), Rodrigues et al. (2017), and Mosser et al. (2013) provided a correction factor to Av (fav), which is calculated 
based on the theoretical Av from fitting radial mode frequencies. This corrects the scaling relations for wide mass 
and metallicity ranges and also for different evolutionary phases, from the ZAMS to the upper RGB. Including these 
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correction factors, the masses and radii reached general agreement for stars in the EB systems (e.g. Brogaard et al. 
2018). More details about the efforts on this topic can be seen in the review by Hekker (2020). 

In this paper, we use masses and radii determined with the detailed modeling to correct the scaling relations. 
Although theoretical models do introduce potential systematic uncertainties, the advantage of using modeling solutions 
is the coverage of a wide range of stellar parameters and evolution stages. For instance, Bellinger (2019) used masses 
and radii previously determined from fits to theoretical stellar models (Silva Aguirre et al. 2015; Davies et al. 2016a; 
Lund et al. 2017a; Silva Aguirre et al. 2017) to calibrate the scaling relations for Kepler dwarfs. We aim to apply this 
method to Kepler red-giant oscillators. We calculate radial model frequencies to determine stellar masses and radii 
and then use inferred values to calibrate the scaling relations. The rest of this paper is organised as follow. Section 
2 summarises the target selection and the measurement procedure of oscillating radial modes; Section 3 discusses our 
theoretical modeling and the fitting method; Section 4 demonstrates results of estimated stellar parameters; we correct 
the scaling relations with model-determined masses and radii in Section 5; and we close with a summary of outcomes 
in Section 6. 


2. DATA REDUCTION 
2.1. Target Selection 


We started with Kepler RGB stars classified by Hon et al. (2017), who used a machine-learning tool to distinguish 
between RGB stars and core-helium-burning stars, based primarily on the gravity-mode spacings (Bedding et al. 2011; 
Mosser et al. 2012b). According to Hon et al. (2017), their classification achieved an accuracy of up to 99%. These stars 
were studied by Yu et al. (2018b), who measured their global seismic parameters (Av and Vmax) using the SYD pipeline 
(Huber et al. 2009b). We also included six oscillating red giants in EB systems, which were previously studied by 
Li et al. (2018), as the benchmark stars in this study. We cross-matched Kepler RGB sample with spectroscopically 
measured atmospheric parameters (T.g and [M/H]) from two large-scale spectroscopic surveys: APOGEE DR16 
(Ahumada et al. 2020) and LAMOST DR5 (Xiang et al. 2019). We only selected stars within a metallicity range of 
[M/H] = -0.3-+0.3 to fall within the parameter range of our model grid (-0.5-+0.5). Note that the [M/H] values are 
not given in LAMOST DR5, and we calculated them from [Fe/H] and [a/Fe] values using the following formula, given 
by Salaris & Cassisi (2005): 

[M/H] = [Fe/H] + log(10/^/*c19.694 + 0.396). (5) 
The cross-match resulted in 4,127 Kepler RGB stars, comprising 2,626 APOGEE and 3,326 LAMOST targets, with 
1,825 common sources. We show the sample on the Teg-vmax diagram in Figure 1. The sample covers the lower to 
the upper RGB, with a Vmax range from 10 to 270 uHz. We also inspected the systematic differences in effective 
temperatures and metallicities between the two surveys. As shown on the right of Figure 1, there is an approximately 
uniform offset in Tyg, which can be described as Tog ApoGEE = Tet,LAMost + 41K. The systematic offset in the [M/H] 
values follows a linear function. We fitted the data and obtained the relation as [M/H]apocrr = 0.20|M/H]rAMosT 
+ 0.023 dex. We show in Section 5 that this systematic offset has an impact on the correction of Vmax scaling relation. 


2.2. Measurement of Radial Oscillation Mode Frequencies 


We first identified the radial modes for the aforementioned stars. The asymptotic relation for p modes can be used 
to predict the mode frequencies on the basis that modes of the same l-degree and sequential n are separated by about 
Av (Tassoul 1980; Gough 1986): 


v(n,l) = Av (r + ! + ) — vol, (6) 


where e defines an offset and dv is the small spacing between modes of different l-degrees. We made use of this 
regularity to determine approximate locations of the radial modes, which means determining e for each star. We 
used the background-corrected power spectra produced by Yu et al. (2018b) and selected data around Vmax in a 
range of 5 times the FWHM (full width at half maximum) of the power excess. We then obtained the so-called 
collapsed power spectrum where the power was added from consecutive slices with width of Av. This collapsed 
power spectrum was cross-correlated with a template comprising three peaks, constructed according to the asymptotic 
equation, corresponding to the ridges of | = 0, 1 and 2 modes. We adopted a value for óvo; that was a fixed fraction 
of Av, set by a fit to real stars (Huber et al. 2010). The value of e was then obtained by maximising the correlation 
coefficient. The locations of radial modes were restricted to be in the region (Vn =0/Av mod 1) € [e — 0.06, e + 0.12]. 
Fig. 2 shows this identification procedure for an example star, KIC 8410637. 
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Figure 1. Left: The Kepler RGB stars on the Tef—Vmax diagram. Grey dots are 50,000 MESA models randomly selected from 
our grid. Blue and red dots represent the 2,626 APOGEE and 3,326 LAMOST tragets. Observed Vmax values are from Yu et al. 
(2018b). Right: Differences in effective temperature (upper) and metallicity (bottom) of the 1,825 common stars. Red dashes 
indicate the systematic offsets, which is constant (Te#,apocEr = Teg,nAMOsT + 41 K) for effective temperature and a linear 
function ([M/H]apocrr = 0.20|M/H]rAwosr + 0.023 dex) for metallicty. 


Having located the radial modes, we fitted Lorentzian profiles to these regions to extract the radial mode frequencies: 


MU) = Bü) it) 3 Lat) = BO): Mm Los UU (7) 


using the fitting method as described by Li et al. (2020b). We used an H1 (odds ratio) approach to evaluate the 
quality of each fit (Appourchaux et al. 2012; Corsaro & De Ridder 2014; Davies et al. 2016b; Lund et al. 2017b). We 
retained those modes with a Bayes factor larger than 3, equivalent to strong detection (Kass & Raftery 1995). As an 
example, we listed model frequencies, uncertainties of KIC 8410637 in Table 1. This star is the red-giant component 
in a classified EB system. Its radial mode frequencies were previously measured by direct fitting (Li et al. 2018). The 
comparison shows that the model frequencies from our pipeline well agree with previous results and we find one extra 
mode at low frequency. We also found good consistencies in extracted modes for five other red-giants components 
studied by Li et al. (2018). 

The number of identified radial modes is 6 on average, reaching up to 10 for low-luminosity RGB stars and decreasing 
with the evolution. For high-luminosity red giants with Vmax below 20~Hz, most stars have only 3-4 modes. We note 
that a few of the stars may have incorrect results since they were returned from an automatic pipeline. They do not 
affect population studies because of their small number. However, to use frequencies for a particular star, they should 
carefully vetted. The source code for the mode measurement and lists of identified frequencies for the all stars are 
available on https:/ /github.com/litanda/SolarlikePeakbagging. 


3. THEORETICAL MODELING 
3.1. Model Grid and Input Physics 


We constructed a stellar model grid that ranges from 0.76 to 2.20M; with a mass step of 0.02Mo. The grid 
includes four independent model inputs: stellar mass (M), initial helium fraction (Yi), initial metallicity ([M/H]), 
and the mixing-length parameter (apr). Ranges and grid steps of these four model inputs are given in Table 2. We 
started computation of each stellar evolutionary track from the Hayashi line and terminated on the upper RGB when 
logg > 1.5dex. Note that we did not include models undergoing core helium burning. 

We used the stellar evolution code MESA (Modules for Experiments in Stellar Astrophysics, version 12115, Paxton 
et al. 2011, 2013, 2015, 2018) and the oscillation code GYRE (version 5.1, Townsend & Teitler 2013) to compute the 
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Figure 2. Measurement of radial modes for KIC 8410637. a) The power spectrum in échelle format. b) The collapsed power 
spectrum. c) The template spectrum. d) The cross-correlation spectrum between b) and c). The red cross shows the maximum 
in the cross-correlation spectrum, which gives the value of e and identifies the ridge of radial modes. We used this to define the 
frequency ranges for the radial modes (shown as black dashed lines) and then fitted the strongest peaks within these ranges. 
'The fitted frequencies are marked by squares. 


model grid. We adopted the solar chemical mixture [(Z/ X)o = 0.0181] provided by Asplund et al. (2009). The initial 
chemical composition was calculated by: 


log(Zinit /Xinit) = log( Z/ X)o + [M/H]ini. (8) 
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Table 1. Extracted radial modes of an example star KIC 10513837. 


v (Hz) cv (uHz) Ref.^ 
37.19 0.02 - 
41.62 0.02 41.63 + 0.03 
46.28 0.02 46.29 + 0.03 
50.85 0.02 50.86 + 0.02 
55.55 0.06 55.53 + 0.05 
60.22 0.05 60.37 + 0.03 


^ Mode frequencies from Li et al. (2018). 


Table 2. Input Ranges and Grid Step of the Model Grid. 


Input Parameter Range Increment 


From To 


M/Mo 0.76 2.00 0.02 
[Fe/H] -0.5 0.5 0.1 
Yinit 0.24 0.32 0.02 
QMLT 1.7 2.5 0.2 


We used the MESA p-T tables based on the 2005 update of the OPAL equation of state tables (Rogers & Nayfonov 
2002) and OPAL opacities supplemented by the low-temperature opacities from Ferguson et al. (2005). The Eddington 
photosphere was used for the set of boundary conditions for modeling the atmosphere. The mixing-length theory of 
convection was applied, where amur = £urr/ Hy is the mixing-length parameter. We considered convective overshoot- 
ing at the core, the H-burning shell, and the envelope. The exponential scheme by Herwig (2000) was applied and the 
diffusion coefficient in the overshoot region given as 


(9) 


2 — 
Dov = Deonv,0 exp (- E) : 


(fo + fov) Hp 


Here, Deonv,o is the diffusion coefficient from the mixing-length theory at a user-defined location near the Schwarzschild 
boundary. The switch from convection to overshooting was set to occur at rg. To consider the step taken inside the 
convective region, (fo + fov)Hp was used, where fo equals 0.5fov. Following Magic et al. (2010), the overshooting 
parameter is mass-dependent following a relation as fo, = (0.13M — 0.098)/9.0, and a fixed fov of 0.018 was adopted 
for models above M = 2.0 Mg. For a smooth convective boundary, we also applied the MESA predictive mixing scheme. 
The mass-loss rate on the red-giant branch with Reimers’ law was set as 7 = 0.2, which is constrained by the seismic 
targets in old open clusters NGC 6791 and NGC 6819 (Miglio et al. 2012). Atomic diffusion was only considered for 
models below 1.1M at the main-sequence phase (it is turned off when central hydrogen fraction goes below 0.01). As 
discussed by Mowlavi et al. (2012), the mass value of 1.1M is a proper switching point, above which the surface helium 
and heavy elements would be rapidly drained out of the envelope if atomic diffusion were applied. The switch of atomic 
diffusion causes a discontinuity in the surface metallicity as well as other global parameters during the main-sequence 
and early subgiant stages, but most of the effects are not significant at the RGB phase due to the redistribution of star 
interiors (Michaud et al. 2010; Dotter et al. 2017). However, element diffusion also changes the main-sequence lifetime 
and hence a discontinuity in star age is expected for the RGB models at 1.1 Mo (e.g. Fig. 11 Mowlavi et al. 2012). 
To investigate the systematic uncertainty, we compared stellar ages at given mean density on RGB for evolutionary 
tracks with and without atomic diffusion. The comparisons shows an 8% difference for 1.1M evolutionary tracks. 
The age difference is strongly mass-dependent, going down to 1% at 1.2M 5 and is approximately zero at 1.3Mo. We 
hence expect a systematic uncertainty of a few per cent in inferred ages in the mass range of 1.1-1.2 Mo. Compared 
with the typical uncertainty (16%), our estimated ages for stars about 1.1Mọ could be biased to some extent for the 
choice of atomic diffusion, but this systematic difference is negligible above 1.2 Mo. 


3.2. Fitting Method 
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We used the Maximum Likelihood Estimate (MLE) method for fitting the models to the observations. In Figure 3, 
we demonstrate our fitting procedure with the example star. We first searched for models in the range 0.8-1.2 times 
Av. Note that the Av scaling relation was used in this step and the +20% range was large enough to cover its 
systematic uncertainty. We fitted the observed T.g and [M/H] to models and calculated the non-seismic likelihood 
(Lnon—seismo) With the following equation: 


(Tee nbs Nn Tef mod) ([M/H]obs = [M/H] moa)? 


2 2 
OT ce. obs 201M Hla, 


Laon—seismo — exp (10) 
where the subscripts ‘obs’ and ‘mod’ stand for the observed and model values. Note that, because our grid resolution 
for [M/H] is 0.1 dex, an uncertainty of +0.1 for [M/H] was used when a star's observed uncertainty is smaller than 
this value. We selected models with Lyon—seismo > 0.0001 for the subsequent seismic analysis. 

We corrected the surface term in oscillation frequencies of each selected model using the two-term formula given by 
Ball & Gizon (2014). While inspecting the fits of frequencies, we noticed that some bad models whose frequencies 
deviated strongly from the observations could be made to fit well after an implausibly large surface correction. To avoid 
this, we set up a threshold for the surface term when selecting models for the seismic fitting. Although the surface 
term of a certain model is unknown before fitting, we can learn from previous studies to roughly estimate its range. 
According to results on Kepler main-sequence and RGB stars (Compton et al. 2018; Li et al. 2018, 2021), relative 
surface corrections at Vmax (OVcorr(Umax)/Ymax) for most stars are in a range from —0.005 to —0.01. Considering the 
potential systematic differences between theoretical models, we used a larger range as the threshold, namely 0 to 0.02. 
After correcting the surface term, we checked the relative surface corrections at Vmax and only kept models fitting this 
threshold criterion. 

We next fitted individual radial mode frequencies and calculated a seismic likelihood of each model for as 


, =i 2 
Leeismo = Į [exp | (Vi obs Vi mod) | . (11) 


(Fh ons F T) 

The subscript i refers to the mode frequencies. As discussed by Li et al. (2020a), the systematic uncertainty in model 
frequencies cannot be neglected in the detailed modeling because it could be comparable to or larger than the observed 
uncertainty. For this reason, the fit was done in two steps. In the first step, we only considered the observed frequency 
uncertainty (0,,555) in the likelihood function and found the best-fitting model. We used the median value of frequency 
offsets between the observations and this best-fitting model as the systematic uncertainty (sys). In the second step, 
we fitted to individual mode frequencies again but included the osys term in the likelihood function. Multiplying the 
non-seismic and seismic likelihoods gave us the final likelihood as L = Lnon-seismo ^ Lseismo- We used the marginal 
likelihood distribution to measure the 16%, 50% and 84% cumulative values to determine global parameters (i.e., mass, 
age, radius, surface gravity, and luminosity) and their uncertainties. 


4. RESULTS 
4.1. Five Red Giants in Eclipsing Binaries 


The accuracy of estimated masses and radii is critical because we use them to calibrate the scaling relations. We 
hence examined the accuracy of modeling solution with five red giants in EB systems, whose masses and radii have 
been previously determined with dynamical modeling. We found APOGEE spectra available for all five red giants 
and LAMOST data for three of them. The comparison is demonstrated in Figure 4. We used the dynamical results 
determined by Frandsen et al. (2013), Gaulme et al. (2016) and Brogaard et al. (2018), and the average value was 
taken when there was more than one estimate for the same star. With the APOGEE spectra, our model-based masses 
and radii agree within lo for three out of the five stars, 1.50 for KIC 9970396, and 2e for KIC 4663623. The reason 
for the large offset on KIC 4663623 could be the very different metallicities adopted in our fitting ([M/H] = 0.08) and 
the dynamical modeling (IM/H] = —0.13). Estimates for the three stars with the LAMOST data are slightly more 
discrepant but still consistent with dynamical determinations. In Table 3, we summarise estimated masses and radii 
given by the detailed modeling, the scaling relations, and dynamical analysis. Compared with the scaling relation, the 
fits to radial mode frequencies significantly improve the accuracy of mass and radius estimates. 


4.2. Estimated Stellar Parameters 


KIC 8410637-APOGEE 
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Top 10 models for KIC 8410637 
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Figure 3. The fitting procedure for the example star (KIC 8410637). Top left: Fitting models to spectroscopic constraints (Teff 
and [M/H]) in the range 0.8-1.2 Av and selecting those with Lnon—seismo > 0.0001 for the subsequent seismic analysis. Grey dots 
represent the all models, colour-coded dots are selected models and colour indicates the non-seismic likelihood (Lon. seismo)- 
Top right: Fitting observed mode frequencies. We demonstrate the fits on the échelle diagram. Filled black squares are observed 
mode frequencies, open blue and red circles represent theoretical mode frequencies before and after the surface correction. The 
échelle diagram is plotted twice to fully present the radial mode pattern, and the blue dashed line indicates the large separation. 
Middle and bottom rows: probability distributions for mass, age, surface gravity, and radius. 
indicates the median of each distribution, and two black dashed lines mark the 16% and the 84% cumulative values (1l-c 
uncertainty). 
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Figure 4. Masses and radii determined by fitting radial modes versus those from dynamical modeling, for five red giants in EB 
systems. The black dashed lines show the 1:1 relation. 
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Table 3. Masses and radii of five red giants in EB systems derived from our modeling, the scaling relations, and dynamical 
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modeling. 
KIC MFit RFit Msn Han Mpyn Rpyn Ref. 
(Mo) (Ro) (Mo) (Ro) (Mo) (Ro) 

4663623 1.567019 9.947932 1.874013 10.78+0.28 1.36+40.09 9.7+0.2 G16 

7037405 1.167008 13.55*037 1.38+0.10 14.7640.41 1.25+0.04 141+0.2 G16 

1.170+ 0.020 14.000 0.093 B18 

8410637 1.544034 10.7073 33 1.78 +0.19 11.4540.56 1.56 0.03 10.7 + 0.1 F13 

9540226 1.324078 12.89*022 1.58+0.11 14.05 +0.37 1.33 0.05 12.8+0.1 G16 

1.378 + 0.038 13.060 + 0.160 B18 

9970396 1.08*008  7.69*015 1.43+0.08 8.68+0.17 1.14 0.03 8.0 + 0.2 G16 

1.178 + 0.015 8.035 +0.074 B18 

References: F13 — Frandsen et al. (2013), G16 — Gaulme et al. (2016), B18 — Brogaard et al. (2018) 
Table 4. Observed constraints and model-inferred stellar parameters (full table available on-line). 
Observed Constraints Estimates 

KIC Source Ter | [M/H] Av Vmax M T R log g 

(K) (dex) (uHz) (uHz)| (Mo) (Gyr) (Ro) (dex) 
8161920 APOGEE 4551 -0.046 3.65 31.9 | 1.144998 7.7133 11.261922  2.393*0007 
5358184 APOGEE 4632 0.104 6.23 62.4 | 1.147903 9.7472 7.877008 2.70170 Opa 
10793872 APOGEE 4780 -0.229 14.86 181.7 | 0.944993 14.6128 4.197008  3.17*000 
8815290 APOGEE 4821 -0.112 446 45.5 | 1.724919 1o9t05  1138'025  2.56*0008 
7749249 APOGEE 4703 0.099 6.66 71.8 |1.22'005 5.1752  'vT2*017  2,751*0007 
8161020 LAMOST 4560 0.240 365 319 | 1.14406% 8.6125 11.251023 2904'0007 
5358184 LAMOST 4573 -0.018 6.23 62.4 | 1.14790? 9.08113 7.872906 —2.702*0 001 
10793872 LAMOST 4925 0.110 14.86 181.7 | 0.945963 102'25. LIGUE 2.171900 
8815290 LAMOST 4749 -0.195 446 45.5 | 1.684013 24755 11.261935 2.55810 000 
7749249 LAMOST 4681 0.117 666 71.8 | 1.241998 4.8152  7.76t917 3.752*0 007 


We used the fitting method described in Section 3.2 to estimate stellar parameters (mass, age, surface gravity, and 
radius) for the 4,127 Kepler RGB stars. We first removed three types of stars from final results: i) stars with three 
or fewer identified radial modes, because the fits will largely affected by the uncertainty in surface correction; ii) stars 
with estimated masses below 0.85 or above 2.10 Mọ, because those values are close to the edge of the grid; iii) stars 
with fewer than 100 models for the seismic analysis, because their observed constraints are either close to the edge or 
outside the parameter ranges of the model grid. We then examined age determinations to identify those targets with 
unrealistically large estimates. These targets may bias the calibration for the scaling relation, given that age correlates 
strongly with mass for RGB stars. We found 61 out of 5,350 estimated ages that are larger than the universe age (13.8 
Gyr) by more than 2-c. We inspected those targets and found some of them do not have reliable Tog measurements 
and others could be core-helium-burning stars according to their power spectra. The fraction (~ 1.1%) of poor age 
determinations generally agrees with the accuracy rate of classifications reported by Hon et al. (2017). We then 
removed the targets with poor age estimates from the sample. We finally obtained 5,289 modeling solutions for 3,642 
stars, including 2,498 APOGEE and 2,791 LAMOST targets with 1,647 common sources. We list observed constraints 
and estimated parameters in Table 4 (the full table is included in the supplement). 

'To inspect the precision of inferred parameters, we show histograms of the uncertainties in Figure 5. We obtained 
similar average precision for the LAMOST and APOGEE samples, but the LAMOST sample spreads over a larger 
range. The typical precision (the median of the distribution) for both samples is 4.596 for mass, 1696 for age, 0.006 
dex for log g, and 1.796 for radius. 

We compared inferred values determined with the APOGEE and the LAMOST spectra to investigate any systematic 
differences caused by the spectroscopy. In Figure 6, we illustrate comparisons for the 1,655 common targets and find 
good consistency at the sample level: the average differences for the four parameter are all around zero, with standard 
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Figure 5. Fractional uncertainties of the estimated parameters for stars with LAMOST and APOGEE spectroscopic constraints. 


deviations that are similar to the typical uncertainties. The large offsets on some stars are due to the very different 
effective temperatures given by the two surveys. The results indicate that our fits are mainly determined by seismic 
mode frequencies, while the effective temperature and metallicity are relatively weak constraints. 


5. THE SCALING RELATIONS FOR RED GIANTS 


Our model-based masses and radii of the Kepler RGB star sample can allow testing of the scaling relations over a 
wide parameter range. We have shown that our estimates for stellar masses and radii are accurate, and hence they 
are good reference values for correcting the scaling relations. We start by inspecting the ratio between the modeling 
determinations and the scaling relations. As illustrated in Figure 7, we find that the scaling relations overestimate 
masses and radii for most of the stars. The density distributions give median factors of ~85% for the mass and ~93% 
for the radius. These values are consistent with those found for EBs (Gaulme et al. 2016; Brogaard et al. 2018; Themefil 
et al. 2018). In Figure 7, we see that the ratio weakly depends on mass and this may correspond to some systematic 
offset caused by the effective temperature or the metallicity. On the other hand, the ratio clearly varies with stellar 
radius, indicating the discrepancy depends on the evolutionary phase. 

We wish to correct the scaling relations by combining the modeling solution and observations. Rather than correcting 
Eq. 3 and 4 directly, we prefer to work with the Vmax and the Av scaling relations because they are more fundamental. 
The scaling relation for Vmax is based on the assumption that Vmax in solar-like oscillators is a fixed fraction of 
Vac(Brown et al. 1991; Kjeldsen & Bedding 1995). The acoustic cutoff frequency (Vac) is the frequency above which 
acoustic waves are not reflected at the surface and it depends on the effective temperature and the surface gravity as 
Vac X JT la Regarding the Av scaling relation, the large frequency separation Av in the asymptotic approximation 
equals the inverse of the sound travel time through the star, which results in Av c y M/R’ (Ulrich 1986). The Av 
value is hence proportional to the square root of the mean density (p) of the star. 


5.1. Correcting the Vmax Scaling Relation 


We calculated the scaling Vmax using model Ter and logg, and used the same fitting method to constrain it. We 
refer to this estimated scaling Vmax aS vmax,sn. Note that we did not use the measured Vmax as an input constraint in 
the fitting, and hence the determination of Vmax,sr is independent of the observed value. 
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Figure 6. Comparisons between estimated parameters with LAMOST and APOGEE spectroscopic constraints. (In the top 
left panel, some data points present stripe-like features. These stars are low-luminous red giants, whose mass loss is negligible, 
and hence their estimated masses are about the same as the initial masses, which are uniformly spaced. Other high-luminous 
stars have signficant mass loss, and hence their inferred masses distribute more randomly.) 


From the comparison between observed Vmax and Vmax sR, as demonstrated in Figure 8a, we find no significant 
systematic offset. However, the scatter (the standard deviation is +2.5%) is much larger the typical observed un- 
certainty, indicating potential dependencies on other stellar parameters or that Vmax uncertainties are substantially 
underestimated. We therefore investigated the correlations with log g, Teg, and /M/H] in Figures 8b, 8c, and 8d. The 
reason to consider the metallicity is that it could impact near-surface properties of solar-like oscillators, as suggested 
by hydrodynamic simulations (e.g. Trampedach et al. 2014). We see that the Vmax offset does not obviously depend 
on logg, but clearly correlates with Tyg and [M/H]. The same trend is seen in both the APOGEE and LAMOST 
samples. From the perspective of physics, the correlations with the effective temperature may be real because the 
thermal structures at the near-surface region are very different between RGB and main-sequence stars, although there 
is no theory that quantifies the difference in Teff-Vmax relations between the two types of stars. However, we suspect 
that this correlation is caused by the systematic offsets in the spectroscopic measurements because the dependency on 
Teg is slightly stronger for the LAMOST data than for the APOGEE data. The same applies for the metallicity, in 
that the correlation is likely a combination of the physical and the instrumental effects and we cannot separate them 
because of the lack of theoretical predictions. We will leave this as an open question for future studies. 

To correct for the systematics discussed above, we added a metallicity term to the Vmax scaling relation. We 
calibrated the revised Vmax scaling relation with model-based log g, observed Vmax, Teg, and [M/H]. We used the SciPv 
CURVE_FIT module to do the fits and obtained corrected Vmax scaling relations as follows: 


—0.459 ME 
la ee ( Ten ) (10/8) ur (APOGEE) (12) 
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Figure 7. Ratios between estimated masses and radii from our modeling and from the scaling relations. Histograms show the 
stars with spectroscopic constraints from LAMOST and APOGEE. 


and 


. —0.421 —0.039 
Vmax = JFit ( Teff ) (10m (LAMOST). (13) 
Vmax,© go Teff,© 


We obtained the same power law for the g term as the original relation, but the exponent of Teg is slightly different 
from the original (—0.5). As discussed above, the difference does not necessarily mean that the RGB stars have a 
different Tef—-Vmax relation from the dwarfs. For instance, when we considered a global temperature shift of —70K 
in the APOGEE sample, as found by Hall et al. (2019), the fit gave a very similar power law to the original scaling 
relation. The additional term of [M/H] is a secondary term in the formulae, which causes a ~4% change to the Vmax 
value when [M/H] varies from -0.3 to 0.3 dex. Note that the LAMOST [M/H] values are based on the [Fe/H] and 
[a/Fe] measurements. The choice of formula may bring some systematic uncertainties and we hence tested the fits 
with another widely-used formula, namely [M/H] = [Fe/H] + [a/Fe]. The [M/H] values from this formula are higher 
than those calculated with Eq. 5 by 0.02 dex, on average, which is significantly smaller than the typical observed 
uncertainty. Fitting the scaling relation with these [M/H] values returns very similar results to Eq. 13. The exponent 
is —0.427 for the Teg term and 0.040 for the [M/H] term. The similarity indicates the choice of [M/H] formula does 
not significantly affect the new scaling relation. 


5.2. The Av Scaling Relation 


We studied the Av scaling relation following the same strategy. We used the mean density of models to calculate the 
large frequency separation from the scaling relation (Eq. 1) and constrained it based on model fitting. We refer to this 
as Avgn. We first examined the systematic difference between observed Av and Avsg. Note that the observed Av 
is from Yu et al. (2018b) rather than from individual radial modes. As shown in Figure 9a, there is a ~4% deviation 
on average, and the offset varies from —6% to —296. We then inspected the dependencies of Av/Avsgr on the mean 
density, effective temperature, and metallicity, as demonstrated in Figures 9b, 9c, and 9d. 

'The ratio clearly correlates with mean density and also with the effective temperature, and slightly depends on 
the metallicity. It is not clear which parameter has the true effect because these three parameters are correlated 
on RGB. As for its physical significance, we note that Av is determined by the mean density rather than by any 
surface properties. As pointed out by Miglio et al. (2012), the differences in internal temperature (hence sound speed) 
distributions could cause a noticeable difference in Av between different type of stars, despite them having the same 
mean density. Their comparison between hydrogen-shell-burning and helium-core-burning stars with similar mean 
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Figure 8. Dependencies of systematic offsets in Vmax (Vmax/Vmax,SR) ON Vmax, Teg, log g,and [Fe/H]. The APOGEE sample is 
presented on the top and the LAMOST sample is shown at the bottom. Filled triangle is the median of binned data and the 
error bar indicates the the standard error of the median, which is calculated as the standard deviation over the number of data 
point in a bin (o/n). 


densities in Kepler open clusters presented a ~3.3% difference in average Av. As a consequence, the difference in 
Av-p scaling relation between RGB and main-sequence stars is expected because of their different internal structures. 
The correlations with effective temperate and metallicity arise because these two surface properties relate to the mean 
density: for stars on RGB, the effective temperature strongly correlates with mass, and hence with mean density, at 
given radius. Moreover, high-resolution spectroscopic surveys have shown a clear T«g-[M/H] relation on the RGB (e.g. 
Fig.2 Tayar et al. 2017). For this reason, we retained the functional form of Av scaling relation and did not include 
the effective temperature or the metallicity in the new scaling relation. 

We calibrated the Av-p scaling relation using observed Av and modeling-inferred mean density (pfit), obtaining 


the following results: 

= 0.507 

Av E (=) , (14) 
Ave Po 


Note that this function is for both the APOGEE and LAMOST samples because there is no systematic differences in 
inferred mean densities between the two at the sample level. Compared with the scaling for Vmax, the scaling for Av 
has a significantly larger systematic offset and the corrections are especially important for RGB stars. The correction 
to the Av scaling relation also matches the evolution-dependent systematic offsets found in previous studies (e.g. 
Sharma et al. 2016). We see that the power law exponent relating the large separation and the mean density for RGB 
stars is larger than that for dwarfs. This causes Av offsets (in fractional terms) to grow exponentially with decreasing 
Av, and hence leads to larger systematic offsets for more evolved stars. We also note that a departure from the 0.5 
exponent in the Av scaling relation has also been seen in 6 Scuti stars (e.g., García Hernández et al. 2015; Bedding 
et al. 2020). 
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Figure 9. Same as Figure 8 but for Av. 


Our correction to the Av scaling relation is similar to previous studies by White et al. (2011) and Sharma et al. 
(2016), who checked whether the Av from fitting the model frequencies matched the value from the scaling the model 
density. This work is an improvement on that earlier work in that our corrections include the surface correction to 
model frequencies. It would be useful to investigate how much the surface term affect the corrections Av scaling 
relation. To do this, we followed the method described by White et al. (2011) and Sharma et al. (2016) to calculate 
the correction factor, fav. In line with their notation, we defined our Av corrections as f 4,, = Av/Avsn. With this 


notation, we can re-write Eq. 14 as 
A = x 
one (15) 


Avo — Po 


Here, f'a, can be taken as an improved version of fa,, with an additional correction for the surface term. Figure 10 
shows the comparison between fA, and f ^4,, of the best-fitting models for all stars. The differences are roughly uniform 
and the average differences are 1.696 and 1.496 for the APOGEE and LAMOST samples, respectively. If we adopt 
a 1.596 difference in the Av correction factor and transfer this difference to the mass and radius determinations, the 
surface term will cause a systematic offset of ~6% for the mass and ~3% for the radius, which are highly significant. 
We hence conclude that the surface term is worthy of attention in the corrections to the Av scaling relation. 

The surface term depends on properties of near-surface layers and we hence inspected its correlations with the mean 
density, effective temperature, and metallicity. As shown in Figure 11, the difference between the two factors decreases 
with Av, implying that the surface term is smaller for more evolved stars. This trend is consistent with 3D convection 
simulations of the surface effect on the RGB (e.g. Trampedach et al. 2017). The effective temperature also has an effect 
on the surface term according to those simulation results, but it is not obvious in our sample because the effective 
temperature range on RGB is fairly narrow. We also notice a rough correlation with the metallicity. This is expected 
because the metallicity can also change the convection properties at near surface Magic et al. (e.g. 2015). 
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Figure 10. The comparison between the two correction factors to the Av scaling relation (fav and f‘',,) for the APOGEE 
sample. Here, f A4, can be taken as the an improved version of fA, due to an additional correction for the surface term. 
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Figure 11. The difference in the two correction factors to the Av scaling relation (fA, - fav) as a function of four stellar 
parameters for the APOGEE sample. Green squares indicates the rolling median of the data. Results for the LAMOST sample 
(not shown) are similar. 


5.3. Corrected Scaling Relations for Red Giants 


With corrected Vmax and Av scaling relations (Eq. 12, 13 and 14), we derived new scaling relations for determining 
stellar masses and radii. For the APOGEE sample they are: 


Mrit _ ( Ymax Bcc cau (soy (16) 
Mo Vmax,© Avo Teft,o , 


—1.972 0.459 
RFit E Vmax Av Tost (omm) an 
Ro Vmax,© Avo Tefo 


and 
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Figure 12. Differences between modeling-inferred results and corrected scaling relations (represented by the subscript ‘CSR’). 
Both density distributions center around zero with a standard deviation of ~ 5% for mass and ~ 2% for radius. 


For the LAMOST sample they are: 


Mp Vmax 3.0 Ap —3.944 Pe 1.263 MTS 0.117 
Mo H (==) (=) (75) (10 ) ’ (18) 
E —1.972 0.421 
RFit _ ( Ymax Av i Tæ NC (gsm yn TO 
Ro Vmax,C Avo Tet, : 


It worth noting that the above relations are only valid for RGB stars within a parameter range of Tog = 4400 to 
5200 K, Vmax = 10 to 270 Hz, and [M/H] = —0.3 to 0.3 dex. We suggest caution when applying them outside these 
parameter ranges. 

As a test, we computed masses and radii with the new scaling relations and compared them with model-determined 
values. We illustrate the distributions of the residuals in Figure 12. The distributions center at zero with a Half- Width- 
at-Half-Maximum of ~5% for the mass and of ~2% for the radius. These are consistent with the typical estimated 
uncertainties, which are 4.596 and 1.796 for mass and radius. 

Because the new scaling relations are based on modeling solutions, some systematic uncertainties are expected. We 
considered reasonable ranges for the helium fraction and the mixing-length parameter, so the systematic effects from 
these two inputs are covered in the estimated uncertainty. However, the choices of other input physics could bias 
the results to some extent. These include the solar composition, the opacity, the atmosphere boundary condition, 
mass-loss rate, and atomic diffusion. The systematic offset, according to comparisons between seven different pipelines 
for Kepler main-sequence stars (Silva Aguirre et al. 2017), is about ~2% for the mass and ~1% for the radius on 
average. Systematic offsets related to input physics in our work could be different because internal structures of RGB 
stars are very different from those of dwarfs. Future studies are required to properly estimate the model systematic 
uncertainties. 

We calibrated scaling relations for the APOGEE and LAMOST samples separately because of the systematic differ- 
ences in their measurements. Given the large number of APOGEE and LAMOST spectra, the new scaling relations 
should be useful for future studies. It is also worth considering whether to apply these new scaling relations to stars 
observed by other spectroscopic surveys. To test this, we applied the APOGEE formulae to the LAMOST data. As 
determined from Figure 1, we uniformly increased the LAMOST Tog by 44K, and calibrated the LAMOST [M/H] 
using the linear function as 0.20[M/H] + 0.023. We used calibrated LAMOST Tog and [M/H] in Eq. 16 and 17 to 
calculate masses and radii. Comparing with model-based determinations, we found good consistency: the average 
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difference was zero with a spread of 01.396 for the mass and +0.4% for the radius. This result indicates that the 
new scaling relations can be applied with spectroscopic data from other sources, provided proper calibrations to the 
effective temperature and the metallicity are implemented. 


6. CONCLUSIONS 


We studied the solar-like oscillations of 3,642 Kepler RGB stars. We measured radial mode frequencies and used 
them to characterise stars with detailed modeling, including correction for the surface effect. We provide model-based 
estimates for mass, age, radius, and surface gravity for the full sample. For five red giants in eclipsing binaries, we 
found good consistency with masses and radii determined from dynamical modeling. We also showed that system- 
atic differences between the APOGEE and the LAMOST measurements of effective temperature and metallicity do 
not significantly affect the modeling solutions at the sample level because the fits are dominated by seismic mode 
frequencies. 

We used our results to study the systematic offsets in the widely-used Vmax and Av scaling relations. We did not 
find a significant systematic offset in the Vmax scaling relation. However, a clear dependency on observed metallicity 
was seen, and we hence added a metallicity term to the scaling relation. We used the observed Vmax, Ter, [M/H], and 
model-inferred log g to calibrate our revised Vmax scaling relation for RGB stars. We noted that the correction to the 
Vmax Scaling relation depends on the selection of spectroscopic measurements. For the Av-density relation, our results 
showed that the scaling relation overestimates the Av value by ~4% on average. We calibrated a new Av-p relation 
for RGB stars based model-inferred p and observed Av, and showed that the surface correction plays a very important 
part in this calibration. 

Combining our results, we derived revised scaling relations for determining masses and radii for RGB stars. The new 
relations can characterise stars with precision of ~5% in mass and ~2% in radius. We suggest that the new scaling 
relations can be applied to stars observed by other spectroscopy surveys, provided any differences in the effective 
temperature and metallicity scales are properly calibrated. 
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